Astronomy & Astrophysics manuscript no. 10579tex 


© ESO 2009 


March 12, 2009 





A lack of close binaries among hot horizontal branch stars in 

globular clusters ^ 

MSOand NGC5986 

O . C. Moni Bidin\ S. Moehler^, G. Piotto^, Y. Momany'^, and A. Recio-Blanco^ 

<^ 

Departamento de Astronomia, Universidad de Concepcion, Casilla 160-C, Concepcion, Chile 
, - European Southern Observatory, Karl-Schwarzschild-Str 2, 85748 Garching, Germany 

^ Dipartimento di Astronomia, Universita di Padova, Vicolo deH'osservatorio 3, 1-35122 Padova, Italy 

INAF-Osservatorio Astronomico di Padova, Vicolo dell'osservatorio 2, 1-35122 Padova, Italy 
' Laboratoire Cassiopee UMR 6202, Universite de Nice Sophia-Antipolis, CNRS, Observatoire de la Cote d'Azur 



Received / Accepted 

ABSTRACT 



C/3 



[ Context. Recent investigations have revealed a surprising lack of close binaries among extreme horizontal branch (EHB) stars in the 

• , globular cluster NGC 6752, at variance with the analogous sdB field stars. Another puzzling result concerns the derived spectroscopic 

i-^ ' masses for some EHB stars. 

Aims. The present paper extends our study of NGC 6752 to M 80 and NGC 5986, to establish whether the unexpected properties of 
EHB stars in NGC 6752 are also present in other clusters. 
<^ . Methods. Twenty-one horizontal branch stars (out of which 5 EHBs) in NGC 5986 and 31 in M 80 (1 1 EHBs) were observed during 

four consecutive nights. We measured radial velocity variations and evaluated statistical and systematic errors. Temperatures, gravities, 
and helium abundances were also measured. 

Results. By means of a statistical analysis of the observed radial velocity variations, we detected one EHB close binary candidate 
per cluster In M 80, the best estimate of the close binary EHB fraction is /= 12%, and even the lowest estimate of the binary fraction 
among field sdB stars can be ruled out within a 90% confidence level. Because of the small observed sample, no strong conclusions 
^ ' can be drawn on the close EHB binary fraction for NGC 5986, although our best estimate is rather low (/=25%). 

I For the discrepancy in spectroscopic derived masses with theoretical models observed in NGC 6752, our analysis of M 80 EHB stars 

, shows a similar trend. For the first time, we report a clear trend in surface helium abundance with temperature, although the trend for 

L~:> • the hottest stars is still unclear 

CNj ' Conclusions. Our results show that the deficiency of close binaries among EHB stars is now confirmed in two, and possibly three, 

globular clusters. This feature is therefore not a peculiarity of NGC 6752. Our analysis also proves that the strangely high spectroscopic 
masses among EHB stars are now confirmed in at least a second cluster Our results confirm that / could be a function of the age of 
the sdB star population, but we find that recent models have some problem reproducing all observations. 



m 
o 

■ Key words, stars: horizontal branch - binaries: close - binaries: spectroscopic - stars: fundamental parameters - globular cluster: 

individual: M 80, NGC 5986 



>< 

H , 1 . Introduction too thin to sustain hydrogen shell burning, and after He exhaus- 

tion in the core they are expected to evolve directly to the white 
Horizontal branch (HB) stars in Galactic globular clusters are ^^^^^^ sequence witho ut ascending the asym ptotic gi- 
old stars of low imtial mass (0.7-0.9 M^) which, after the ant branch (AGB manque stars. lCkeggio & Renzinil ll9901. EHB 
exhaustion of hydrogen in the stellar core and the ascension ^^^^ ^^^^ ^^^^ extensively observed and also studied in the 
along the red giant branch, eventually ignited helium core Galac tic field, ide ntified as the so-c alled su bdwarf B-type (sdB) 
burning ( [Hoyle & Schwarzschild| [19551 iFaulknd ll966|). The ^^^.^ (iGreensteinll 19711 ^MM^ Meb^\m6 ). although this 
most puzzhng feature of these stars is surely the large va- ^ spectroscopic classification without direct link to the stel- 
riety of HB mophologies in globular clusters (see for ex- evolutionary stage. Given the intrinsic faintness of these ob- 
ample Piotto et al. | |2002^, which is only pardy explained by -^^^^^ ^^^^ ^^^jj spectroscopically poorly observed in globular 
difi-erences in metallicity (thejo-called "second parameter clusters, and many open problem s lack full comprehension (see 
problem", [ Sandage & Wildeyl ll967t IvaiLden Bergh| ll967|). In |Moni Bidin et alllSbt IC^tei^l200l for recent reviews), 
this context, the foremost problem is the presence of ex- 
treme horizontal branch (EHB) stars at the faint hotter end of One of the most evident features of HB stars is the on- 
HBs (Teff >20000 K), eve n in high metalli city clusters Hke set of atmospheric diffusion for temperatures T^ff >11000- 
NGC6388 and NGC 6441 dRich et all ri997h . EHB stars are 12 00 K. This catises ch anges in their photometric proper- 
identified as hot He-core burning stars with an external envelope ties dGrundahl et all Il999h . deficiency of helium because of 

gravitational settling, and strong (solar to super-solar levels) 

* Based on observations with the ESO Very Large Telescope at enrichment of heavy metals do iaspey et al] 1 19851 NGC6397; 

Paranal Observatory, Chile (proposal ID 69.D-0682) iGlaspev etalll 19891: iMoehler etal..,200a NGC6752; iBehr et all 
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1991 iMoehleretalJ llool M13; iBehretaP llOOOl M15; 



Fabbian et all 1200 51. NGC1904; 'Pace et alj |200l, NGC2808) 



Precise calculations confirmed that diffusion should be at work 
in the atmos pheres of these star s and can account for observed 
anomalies dMichaud et al.lll983h . lMichaud et all (l200 8l) recently 
confirmed the role of atomic diffusion in the observed abundance 
anomalies using new stellar evolution models. Spectroscopically 
determined surface gravities are systematically lower than pre - 
dictions for stars in this temperature range (Mo ehleret al.ll995h . 
a problem only partially explai r ied by abundance anomalies 
(lMoehledl200TD . IVink & Cassisil (|2002|) point out that neglect- 
ing the presence of stellar wind can cause measured surface 
gravities to be erroneously lo w, and an enhanced stell ar wind is 
actually the explanation that iMoni Bidin et al. I (120071 hereafter 
M07) proposed for some bright stars showing erroneously low 
masses. It is worth noting, however, that optical and UV obser- 
vations do not support high mass loss rates for field EHB stars 
in general, with the exception of few re lat ively luminous objects 
(iMaxted et al.ll200ll: iLisker et al.ll2005h . iMomany et aP (l2002l 
|2004 >) foundthat at temperatures hotter than 23 000 K HB stars 
deviate again from canonical tracks in the color-magnitude di- 
agram, and they proposed a new onset of diffusion as expla- 
nation of this feature. T he low helium abundances found by 
iMoehler et al.l (l2000l) and lM07l on some hot stars seem to con- 
firm this hypothesis, but the pattern of abundance with tempera- 
ture is unclear 

In canonical models, the mass of the He-burning core is ap- 
proximatively the same for all HB stars, equal to the minimum 
requir ed for core helium flash O.SM©, Schwarzschild & Harm 
1 19621) . but the envelope mass decreases for higher temperature. 
The extremely hot EHB sta rs retain just a very thin inert hy- 
drogen envelope (<0.02Mo, lHeber|[T986h . and must have suf- 
fered an extreme mass loss during their evolution. Many single- 
star evolutionary channels have been invoked to explain EHB 
star formation i n globular c lusters, in cluding interactions with 
a close planet jSokeJ I1998L see also 'Silvotti et al.l 12007'), He 
mixing driven b y internal rotation {Sweigart & Mengel 1979; 
ISweigarj Il997h or by stellar encounters (Sudaetal. 2007), 
dredg e-up induced b y H-shell instabiliti es (von Rudloff et al, 
11988 , but see also De nissenkov & VandenB erg 2003D, close en- 
count ers with a central, intermediate-mass black hole (i Miocchil 
I2007h . and a sub-population of s tars with high helium abun- 
dance (e.g., iD'Antona et al.l 2005h. T he discovery of multiple 
main sequences in c^i Cen (iBedin et al. 2004) and in NGC2808 
dPiotto et al.l l2007h reinforced the idea that in some clus- 
ters there might be a fraction of stars super-He rich, up to 
Y~0.40 ( Norris 200^ iPiotto et al.ll200l ID'Antona et aLll200l 
iLee et all 1200 5). Nevertheless binary models, in which sdB 
stars form through dynamical interactions within binary sys- 
tems, have been very s ucces sful in reproducing observations 
dHanet al.ll2002l l2003l 12007*). and are actually the most pre- 
ferred scenario for field sdB star formation. Indeed, many sur- 
veys hav e shown the existence of a large population of sdB 
binar ie s (iFergu son et al.lll984l: l Allard et al.ll 19941: lUUa & Thejlll 
19981: lAznar Cuadrado & Jeffer vl 120011: iMaxted et all 120011: 



Williams et all i2001c iReed & St iening '200?, 'Napiwotz ki et all 
2004 . Among them, clos e system s with periods shorter than 



10 days play a major role (iMoran et a l. 1999; S affer et al .1119981: 
iHeber et al.i2002l: iMorales-Rueda et al...2003.) . The close binary 
fraction among field sdB stars is certainly high b ut still ill- 
determined, ranging from 70% (I Maxted et "511200 lb to 40-45% 
(Napiwotzki et al. 2004). In this context, it came as a great sur- 
prise that first surveys in globular clusters r evealed a lack of 
close binary systems among the EHB stars dMoni Bidin et al.l 



I2006bi hereafter Paper I). Recently iMoni Bidin etall ( l2008ah 
showed that the best estimate of the close binary fraction among 
EHB stars in NGC 6752 is only f-4%. On the basis of theoreti- 
cal and observational results available in the h terature, they sug- 
gested the presence of a /-age relation. iHanl (12008) supported 
this hypothesis with detailed theoretical calculations, showing 
that the binary scenario naturally predicts a steep decrease of 
close binary fraction with increasing age of the sdB population. 
This seems a further success of Han's models, but the general 
lack of observational data in globular clusters and the uncertain- 
ties on the predicted / values (due to uncertainties on model 
parameters) still requires caution. 

Observations of globular clusters other than NGC 6752 are 
needed to verify that the low fraction of close EHB systems is 
not a peculiar feature of this cluster. Also other types of binary 
systems must be included in future searches. In fact, bo th wide 
binaries dReed & Stiening|2004l;lMorales-Rtieda et al . 2006) and 
system s with very low-mass secondaries (iMenzies & Maran^ 
[l986h are known to exist among field sdBs (although they are 
a minor population there) and it has been shown that they can 
provide an additional channel for the formation of sdB stars. 

In this paper we present the results of a binary search in two 
additional globular clusters, NGC 5986 and M 80. A f irst analy- 
sis of data was presented in iMoni Bidin et al] d2006ah . Here we 
refine that preliminary overview with the correction of system- 
atic effects and a detailed error analysis, and we extensively use 
statistical calculations to better clarify the significance of our re- 
sults. We also present results about atmospheric parameters and 
masses for our target stars. 



2. Observations and data reduction 

We selected 21 HB stars in NGC 5986 and 31 in M80, span- 
ning a wide range in effective temperature. Targets in M 80 were 
divided in two fields, named M80a (17 stars) and M80b (15 
stars), with different slit configurations for multi-object spec- 
troscopy. One star in M 80 (#14327) was accidentally observed 
in both fields (as star la and 12b respectively). The positions 
of the observed stars along the HB of their parent cluster are 
shown in Fig. [ T] whi le ast rometric and photome tric data from 
iMomanv etall d2003l) and iMomanv etal] d2004 are presented 
in Table 12] 

The spectra were collected during four nights of observations 
(June 11th to 14th, 2002) at the VLT-UT4 telescope equipped 
with the spectrograph FORS2 in MXU mode. We employed the 
grism 1400V with 0'.'5- wide slits, and the resulting resolution 
was 1.2 A. The 2400 s exposures were always acquired in pairs 
and subsequently summed, with the exception of one single ad- 
ditional spectrum of NGC 5986 targets during the second night. 
We finally decided to exclude it from analysis because of its too 
low S/N, which caused unreliable measurements. The bias, flat, 
and lamp exposures were acquired during daytime. Just before 
each pair of exposures, a slit image (without grism) was taken, 
which was used to correct the spectroscopic data for instrumen- 
tal effects. The spectral range varied from star to star because of 
different positions of the slit in the mask, but the Hyj line was 
always inside the spectral range. 

During each night, we successfully collected at least one pair 
of medium-resolution spectra for NGC 5986, and the resulting 
temporal sampling was very good. Unfortunately, because of 
strong winds from north, observations of M 80 could not be car- 
ried out during the entire first night and partially the second one. 
The starting time of the pairs of exposures is given in Table [1] 
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Fig. 1. Color-magnitude diagram of NGC 5986 {upper panel) 
and M 80 (lower panel), with targets highlighted. Stars for which 
we derive anomalously hi gh masses are plotted a s full triangles. 
Pho tometric data a r e fromlMomany etal] (l2003h for NGC 5986 
and lMomanv et all (|2004 for M 80. 



During the same run we took two 1350 s exposures of each 
target with grism 600B, for a resulting resolution of 3A, to mea- 
sure atmospheric parameters. Spectra were trimmed at 3600 A 
on the blue side, because of the lack of instrumental response 
and atmospheric transmission. All Balmer lines from to HIO 
were always present in these spectra. 



Field 


nii 

12 13 


ihl 

14 


15 


NGC 5986 


1:18 23:19 
3:23 


3:00 


3:00 


M80a 


6:20 


23:27 
0:53 


4:30 


M80b 




4:31 


23:32 
0:57 



Data reduction was performed with standard MIDA^proce- 
dures. All slitlets were trimmed from the multi-object frames and 
reduced independently. The wavelength calibration (wlc) was 
performed using He and HgCd lamp exposures, fitting a 3'^'' order 
polynomial to the dispersion relation for both grisms. The mean 
rms of this fit for the medium-resolution spectra was 2.46- lO^^A. 
All two-dimensional spectra were corrected for curvature along 
the spatial axi s tracing them with a speci fic MIDAS routine, 
as described in iMoehler & SweigarB (12006), and then extracted 
with an optimum extraction algorithm (Horne 1986). Sometimes 
(usually for brighter stars) this procedure failed, producing noisy 
spectra with irregular continuum, and we opted for a simple sum 
algorithm in these cases. 1400V spectra were rebinned to con- 
stant steps of 0.25 A/pix, and continuum-normalized. Most spec- 
tra of M 80 stars showed a tiny but clear interstellar emission line 
in the core of H^. During data reduction we gave particular at- 
tention to the proper removal of this feature, which could spoil 
the RV variation measurements. 600B spectra were rebinned 
to a larger constant step (0.4 A/pix), corrected for atmospheric 
extinction with the coefficients for La Silla observatory (iTiigl 
and flux-calibrated. The response curve was obtained sep- 
arately for each night, through observations of standa r d star s 
EG274 and LTT3218 with the flux table of Ha muvetai] ( [T994h . 
Finally, on 600B spectra we fitted a Gaussian profile to the core 
of all Balmer lines from to Hg, excluding due to blending 
with the Can H line, and we used the resulting average radial 
velocity to shift the observed spectra to laboratory wavelengths. 



3. Measurements 

3.1. Atmospheric parameters arid masses 

To derive atmospheric parameters by means of Balmer and he- 
lium lines fitting, it is important to know whether the stellar 
atmosphere is affected by diffusion processes. In fact, the pro- 
file of lines under study c an be influenced by helium deple- 
tion and metal enrichment. IMoehler et all (1 19991 l2000l l2003h 
showed that strong Fe ii lines in the region 4450-4600 A are de- 
tectable even at low resolution when diffusion starts up at about 
12 000 K (see for example Fig. 3 in iM07 ). Therefore, spectra 
showing evidence of iron lines or of being hotter than 14 000 K 
(as deduced from their position in the color-magnitude diagram) 
were fit with metal-rich ([M/H] -1-0.5) model spectra, whereas we 
adopted metal-poor models ([M/H] - -1.5) for all other stars. 
In few cases we relied also on higher resolution 1400V spectra, 
summed altogether for a higher S/N, in search for evidence of 
atmospheric diffusion. We kept the helium abundance fixed at 



' ESO-MIDAS is the acronym for the European Southern 
Observatory Munich Image Data Analysis System which is de- 
veloped and maintained by the European Southern Observatory 
^http://vywvy.eso.org/projects/esomidas/j 



4 



C. Moni Bidin et al.: A lack of close binaries among hot horizontal branch stars in globular clusters 



Table 2. Photometric data of program stars. Columns 1: slit 
number Columns 2-6: ID s, coordinates and photometric data 
from ,Momanv et al., (.20031) and lMomanv et al. (.2004.) . 



slit 


ID 


RA (J2000) 


DEC (J2000) 


V 


(U-V) 






hh:mm:ss 








NGC5986 


1 


17512 


15:46:06.215 


-37:50:49.84 


18.320 


-0.052 


2 


17604 


15:46:03.177 


-37:50:41.41 


16.892 


0.756 


3 


17691 


15:46:08.614 


-37:50:29.56 


19.812 


-0.303 


4 


18077 


15:46:05.939 


-37:49:59.20 


17.016 


0.691 


5 


18240 


15:46:07.543 


-37:49:47.01 


16.851 


0.749 


6 


2103 


15:46:10.144 


-37:48:31.16 


17.559 


0.190 


7 


3192 


15:46:11.382 


-37:48:06.34 


18.067 


-0.077 


8 


3560 


15:46:12.748 


-37:47:58.32 


18.434 


-0.099 


9 


4175 


15:46:16.056 


-37:47:44.62 


19.417 


-0.279 


10 


4930 


15:46:18.214 


-37:47:29.39 


17.143 


0.453 


11 


5558 


15:46:11.971 


-37:47:21.61 


17.349 


0.252 


12 


6102 


15:46:08.134 


-37:47:13.91 


17.725 


0.147 


13 


7008 


15:46:13.634 


-37:46:54.66 


16.879 


0.604 


14 


7430 


15:46:11.091 


-37:46:48.42 


17.004 


0.514 


15 


8131 


15:46:11.505 


-37:46:34.93 


19.370 


-0.345 


16 


9049 


15:46:04.640 


-37:46:21.30 


18.292 


-0.105 


17 


9250 


15:46:17.684 


-37:46:09.79 


16.929 


0.673 


18 


10390 


15:46:09.599 


-37:45:51.45 


16.681 


0.554 


19 


11215 


15:46:11.234 


-37:45:30.85 


17.182 


0.446 


20 


11571 


15:46:16.004 


-37:45:17.74 


20.078 


-0.402 


21 


12099 


15:46:09.689 


-37:45:05.36 


19.282 


-0.189 


M80 


la 


14327 


16: 


17:11.476 


-22:59:23.30 


19.086 


-0.386 


2a 


14786 


16: 


17:10.954 


-22:59:06.21 


16.304 


0.596 


3a 


14985 


16: 


17:08.723 


-22:58:59.74 


18.290 


-0.223 


4a 


15200 


16: 


17:06.303 


-22:58:53.05 


18.152 


0.084 


5a 


16389 


16: 


17:12.046 


-22:58:05.56 


18.549 


-0.379 


6a 


16163 


16: 


17:05.336 


-22:58:18.14 


17.724 


-0.117 


7a 


17173 


16: 


17:10.234 


-22:57:37.98 


16.244 


0.744 


8a 


17114 


16: 


17:07.499 


-22:57:42.36 


17.515 


-0.045 


9a 


16707 


16: 


16:56.026 


-22:58:04.77 


20.023 


-0.629 


10a 


17737 


16: 


17:03.158 


-22:57:21.75 


19.359 


-0.518 


11a 


18516 


16: 


17:07.675 


-22:56:45.43 


18.635 


-0.344 


12a 


18110 


16: 


17:00.961 


-22:57:07.79 


17.053 


0.027 


13a 


18992 


16: 


17:08.007 


-22:56:18.24 


16.391 


0.469 


14a 


19040 


16: 


17:06.095 


-22:56:16.68 


18.574 


-0.337 


15a 


18391 


16: 


16:54.498 


-22:56:59.78 


15.842 


0.691 


16a 


19191 


16: 


17:03.050 


-22:56:07.66 


19.124 


-0.483 


17a 


19127 


16: 


16:55.754 


-22:56:16.88 


19.450 


-0.552 


lb 


12304 


16: 


16:58.233 


-23:01:19.25 


19.590 


-0.641 


2b 


14201 


16: 


16:57.762 


-22:59:37.14 


17.898 


-0.171 


3b 


13787 


16: 


16:59.111 


-22:59:53.45 


18.261 


-0.205 


4b 


12663 


16: 


17:01.896 


-23:00:48.22 


19.133 


-0.59 


5b 


12767 


16: 


17:03.678 


-23:00:40.16 


19.712 


-0.625 


6b 


13839 


16: 


17:03.211 


-22:59:48.33 


16.342 


0.596 


7b 


14022 


16: 


17:05.686 


-22:59:38.95 


18.909 


-0.245 


8b 


14387 


16: 


17:06.469 


-22:59:24.31 


17.697 


-0.050 


9b 


12526 


16: 


17:10.297 


-23:00:52.02 


19.520 


-0.688 


10b 


13179 


16: 


17:10.249 


-23:00:12.68 


17.348 


-0.284 


lib 


15682 


16: 


17:08.241 


-22:58:34.56 


16.995 


0.005 


12b 


14327 


16: 


17:11.476 


-22:59:23.30 


19.086 


-0.386 


13b 


15183 


16: 


17:12.947 


-22:58:49.48 


15.866 


0.692 


14b 


14813 


16: 


17:14.702 


-22:59:02.51 


16.427 


0.485 


15b 


15470 


16: 


17:15.921 


-22:58:37.19 


16.034 


0.725 



solar value (log ^ = -1.00) for cool stars (Teff <11000K), 
as the helium lines in their spectra are rather weak. During the 
fitting we verified that helium lines predicted for these targets 
agreed with the observed ones. 



W e computed model atmospheres using ATLAS9 (iKurucd 
1993") and used Lemke's versiorQ of the LINFOR program (de- 
veloped originally by Holweger, Steffen, and Steenbock at Kiel 
University) to compute a grid of theoretical spectra that include 
the Balmer lines H„ to H22, Hei (4026 A, 4388 A, 4471 A, 
4921 A), and Hen Hnes (4542 A and 4686 A). The grid cov- 
ered the range 7 000 K < Teif < 35 000 K, 2.5 <\ogg < 6.0, 
-3.0 < log^ < -1.0, at metallicities of [M/H] = -1.5 and 
+0.5. To establish the be st fit to the observed spec tra, we used 
the routines developed by Beraeron et al. ( 1992) and lSaffer et alJ 
( 1994 ). as modified by Napiwotzki et al. ( 1999), which employ 
a test. The cr necessary for the calculation of is estimated 
from the noise in the continuum regions of the spectra. The fit- 
ting program normalizes model and observed spectra using the 
same points for the continuum definition. was excluded to 
avoid the blended Can H line. The errors in eac h fitting proce- 
dure w ere derived from the;if^ of the fit itself (see lMoehler et al.l 
'1999' for more details), under the assumption that the only er- 
ror source is the statistical noise. However, Napiwotzki (priv. 
comm.) noted that the routine underestimates this statistical er- 
ror by a factor of 2-4. In addition, errors in the normalization 
of the spectrum, imperfections of flat field/sky background cor- 
rection, etc. may produce systematic errors, which are not well 
represented by the error obtained from the fit routine. 

For each star, we measured atmospheric parameters in the 
two 600B spectra separately, and the final results are the 
weighted mean of the two measurements. Errors were multiplied 
by V3 because the fitting procedure assumes each pixel as inde- 
pendent of the others, but when rebinning we oversampled the 
spectra by a factor of three with respect to the dispersion. Stars 
la and 12b in M80, which are actually the same object, were 
studied as if they were two different targets, so we could compare 
the results as an indication of their quality. Results are in excel- 
lent agreement (see Table HJi, despite the target being among the 
faintest ones, showing that the reported internal errors are proba- 
bly realistic. The parameters from the two distinct measurements 
are so similar that we found no need to weight them altogether, 
hence in our analysis we will simply omit star 12b in order to 
assure statistical uniformity to our sample. Here, the choice of 
which star to exclude is quite irrelevant, and we opted to keep 
star la for continuity with RV variation analysis, where the two 
slitlets are not equivalent because of the better time coverage of 
field M 80a (see discussion in 33.2b . 

Masses were calculated from the previously measured atmo- 
spheric parameters, through the equation: 

, M g T M^+BC- 4.14 

log— = log 4 -log- — , (1) 

Mq gQ Tq 2.5 

obtained from basic relations. We adopted the standard values 
Tq=5777 K and loggQ=4.4377. The bolometric correction was 
derived frorneffective temperature through the empirical calibra- 
tion of iFlowerl (II99 6). We adopted an apparent distance modu- 
lus (m--M)v^ l5.96 for NGC 5986 and (m-M)v=15.56 for M 80 
( lHarris|[T996l February 2003 Web version). Errors on mass es- 
timates were derived from propagation of errors, assuming an 
uncertainty of 0. 1 on photometric quantities (distance modulus, 
magnitude and BC). 

Our results are presented in Table |4] Despite the fact that in- 
struments, observing nights, and measurement procedures were 
the same as in M07, the errors are much higher here. This is eas- 
ily explained by the faintness of the targets in these two clusters 
with respect to NGC6752. 



For a description see |http://a400.sternwarte.uni-erlangen.de/~ai26/linfit/lirrfor.htn 
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Fig. 2. Upper panel, absolute value of measured RV variations 
of all stars plotted against the spectral S/N. Lower panel: same 
plot but as a function of stellar effective temperature. 



3.2. Radial velocity variations 

Radial velocity (RV) variations were measured by mea ns of the 
cross-correlation (CC) technique (iTonrv & Davi s 1979). The tar- 
gets are fainter by 1-2 magnitudes with respect to Pap er Ij, and 
spectra were much noisier despite the longer exposure times. 
As a consequence, the measurement procedure presents some 
small difference with respect to previous work, and final errors 
are higher, especially for hot stars. 

We cross-correlated each spectrum with all the others of the 
same star using the IRAF0 task /xcor. Thus, we measured 10 
RV variations for each target in NGC 5986 (5 spectra per star), 
6 in M 80a (4 spectra), and only 3 in M 80b (3 spectra). To de- 
fine the sign of the variation, unnecessary for our goals, the first 
spectrum in temporal order was always assumed as template. 
We cross-correlated the single spectra of each pair before sum- 
ming them, to check that the sum was safe and no RV variation 
occurred in between exposures. Unfortunately, single spectra of 
very hot stars were often of too low S/N for a reliable CC, and 
they were summed bona fide. 

RV variation measurements focused on Hp line with full 
wings. Based on our experience, the CC of such a wide line fails 
to give the correct shift of the spectra if restricted to line cores, 
and results tend to lower values if the adopted interval is too nar- 
row. On the other hand, on larger intervals, the results converge 
to a fixed value, but the line wings provide progressively less 
information and more noise. We performed tests on artificially 
shifted spectra, looking for the best compromise, and eventually 
we adopted the interval 4830-4890 A in our CC. The CC func- 
tion roughly resembles the line profile, and in our case it took 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 




canonical 
polluted (Y-0.33) 



20000 



10000 8000 



Te„ (K) 



Fig. 3. Teff-log(g) plot of observed stars, (full points: NGC 5986, 
open points: M 80). Errorbars on single stars are omitted for clar- 
ity, but the eiTorbar drawn in lower right corresponds to typ- 
ical values ^^f^=0.03, cr(log(g) )^0.07 dex. Small points in- 
dicate results on NGC 6752 from lM07i Zero-age (ZAHB) and 
terminal-age (TAHB) horizontal branch theoretical tracks are 
also indicated, for [M/H]=-1.5 an d canonical models ( Y-0.24) 
and polluted ones (Y=0.33) (see iMoehler et all 120031 for de- 
tails). 



the shape of a Gaussian core with wide wings. We performed a 
Gaussian fit of the central peak to determine its center. For hot- 
ter stars, the l ow S/N (s ee Fig. |2]i required the application of a 
Fourier filter i Brault & White 1971). The resulting CC function 
lost its strong irregularities and was easier to fit, but we always 
verified that RV variation after filtering did not differ from what 
could be (sometimes with difficulties) deduced fitting the unfil- 
tered CC function 

We always tried to confirm the measured RV variation cross- 
correlating weaker spectral lines. On hotter stars this proved 
fruitless, because too low S/N caused unreliable (often even im- 
possible) measurements. Results agreed with H^, but within so 
large eiTors that the confirmation was useless. On cooler stars we 
obtained very good CCs but they just confirmed measurements 
on Wp, without additional information. Hence, we will not ana- 
lyze results obtained with lines other than Hp. 

3.2.1. Corrections on radial velocity variations 

R V variati ons were corrected as in lPaper j The reader is referred 
to lPaper j for an extensive discussion. In brief, we first used the 
[O i] 5577 A sky emission line as a zero-point to coiTect the spec- 
tral shift with respect to the arc lamp observations. The position 
of the forbidden line was determined with a Gaussian fit. After 
this, we corrected RV variations caused by different positions 
of the stars within the slits, using the slit images secured before 
each pair of exposures. 
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Table 3. "Extraction and fit" errors (in km s ' ) for each field and 
different S/N range. 



S/N 




field 






NGC 5986 


M80a 


M80b 


>80 


1.2 


1.2 


1.0 


40-80 


2.1 


3.0 


3.0 


<40 


6.3 


5.0 


5.4 



The corrections to be applied clearly correlate wit h the spa - 
tial Y coordinates on the CCD (see Fig. 5 and 7 of iPaperl) . 
In both steps we preferred to obtain the final corrections from 
a least-square fit as a function of Y, to reduce the additional 
noise added to the final results. The correction procedure on 
field M 80b was straightforward. On NGC5986 and M 80a we 
sometimes found mismatches between the derived corrections 
and RV variations to be corrected, in a way very similar to Fig. 7 
of|Paper I (middle panel). This indicates movements of the mask 
inside its fra me betw een the slit and science images, as akeady 
discussed in Paper J . On the other hand, we found negligible in- 
dication of rotation of the masks (i.e. different slope between 
corrections and RV variations when plotted against spatial coor- 
dinate). It was recently founcQ that the central wavelength for 
FORS2 data varies with the temperature at the telescope fo- 
cus. This fact also could explain the observed shifts of spectra 
on CCD, but we found no clear coiTelation between them and 
focus temperature differences as obtained from frame headers. 
Moreover, no shift was observed a mong sp ectra in field M 80a, 
as among many frames studied by iPaper iL despite the fact that 
temperature changes occurred. We must conclude that in some 
cases we observe small rigid shifts of spectra on CCD not caused 
by different centering of stars in the slits, that can be coiTected 
because th ey affect simultaneously all the spectra in the same 
frame (see iPaper ll for details). They cannot be explained by 
variations of the temperature at telescope focus alone, and move- 
ments of the mask within its frame must play a role, although the 
two causes could act together and sum their effects. 

3.2.2. Errors on radial velocity variations 

The analysis of RV variations in search of binary systems re- 
quires an accurate error analysis, because the crucial point is to 
tell if variations are consistent with random measurement eiTors 
or are an indication of binarity. As discussed later (©, varying 
estimated errors by 10% can noticeably change the probability 
of the datum, up to a factor of two. 

We estimated the error associated to each RV variation as 
the quadratic sum of all relevant sources. The wavelength cali- 
bration error was deduced from the rms of the calibration pro- 
cedure, and it was directly calculated by the calibration routine. 
Its exact value changed from spectrum to spectrum, but the dif- 
ferences were negligible in the final quadratic sum, therefore we 
kept this contribution fixed to the mean value of 1 .5 km s ' . This 
is in good agreement with the eiTor estimated in Paper I for the 
same instrumentation, obtained analyzing calibrated lamp im- 
ages. The error introduced by the correct ion of systematic ef- 
fects 03.2.1t were estimated as in lPaper j, from the uncertainty 
in the sky line peak position, and the scatter of residuals with 
respect to the least-square solutions used for the co rrections . We 
did not find any relevant difference with respect to iPaper II and 

see http://www.eso.Org/observing/dfo/quality/F ORS2/reports/l 
HEALTH/trendj-eport_LSS.lambda_c.T_HC.html 



we adopted a fixed cTsjj^y^l.S km s for the first correction un- 
certainty and cTcon-^-i km s"' for the second one. Wavelength 
calibration and sky-line correction errors were considered twice 
in the quadratic sum, because in each RV variation two spectra 
were involved. 

Additional uncertainties were introduced by the choice of the 
parameters determining the CC function fit and of the spectrum 
extraction, because different extractions caused noise-induced 
differences in the line profile. They were estimated re-extracting 
spectra in slightly different ways and re-fitting CC functions for 
all hot stars and a sample of cooler ones. The dispersion of the 
differences between the measurements was assumed as an es- 
timate of the error. This parameter is very sensitive to spectral 
noise, and we grouped stars in each cluster in three ranges of 
S/N. The resulting errors are given in Table[3] In Fig.|2]we plot 
the absolute value of RV variations for all stars, as a function 
of their S/N and temperature, excluding the three binary candi- 
dates discussed in ^(two EHB and a cool target). The trend of 
decreasing dispersion with increasing S/N is clear. 

3.3. Absolute radial velocities 

Absolute RVs were measured on 1400V spectra by means of 
CC with synthetic spectra of appropriate temperat ure and grav- 
ity drawn from the library of iMunari et al.l (l2005h . We verified 
that the template metallicity had negligible effects on the results. 
Each measurement was corrected to make the [O i] 5577 A sky 
line coincide with its laboratory wavelength, then the weighted 
mean was computed to derive the final absolute RV. The error re- 
sulting from the weighting procedure was unrealistically small, 
and for the final error we adopted the dispersion of the single 
measurements. Results are shown in Table |4] 

Within errors, all stars show an absolute RV compatible with 
parent cluster. We can conclude that NGC 5986 targets are most 
likely cluster members. For M 80 stars we can just state that RV 
does not disprove cluster membership. In fact, because of the 
very low cluster RV (8.2 km s"', Harris 1996, February 2003 
Web version), we cannot distinguish members from foreground 
Galactic disk stars, which are expected t o contaminat e the field 
at such low Galactic latitudes (b - 19!5. lHaiTisll 19961 February 
2003 Web version). Indeed, we have reasons to suspect that some 
cool stars are main sequence foreground objects (as discussed 
in ®, but their absolute RVs do not noticeably differ from the 
cluster value. The RV of star #14327 in M 80 was measured sep- 
arately in both fields, but results are identical in value and error. 

4. Results on atmospheric parameters and masses 

The 600B spectrum of star #16707 in M 80 fell between the two 
FORS2 chips, and we could not measure its atmospheric param- 
eters. Its 1400V spectra were well inside the second chip, be- 
cause of the offset between the two grisms in the spatial direc- 
tion, and RV variation measurements proceeded normally. We 
deduced its effective temperature from its color, for the only pur- 
pose of its classification in our search for binarity (® . 

Our results on effective temperatures and gravities are plot- 
ted in Fig. |3] where we compare the position of our tar- 
get stars with theoretical models. The zero-age (ZAHB) and 
termi nal-age (TAHB) HE theoretical tracks for [M/H]=-1.5 
from iMoehler et all (l2003h are indicated, both for canonical 
models of normal He-content (Y=0.24) and polluted ones with 
enhanced helium abundance (Y=0.33). The ZAHB and TAHB 
define the region where models spend 99% of their HB lifetime. 
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Fig. 4. Calculated masses of program stars as function of effec- 
tive temperature. Symbols are as in Fig. [3] M80 stars which 
show anomalous masses are plotted as full triangles. Errorbars 
on temperatures are o mitted for clarity. Also the theoretical HB 
dMoehler et al.ll2002l) is indicated. 
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Fig. 5. Measured surface helium abundance for program stars. 
Symbols are as in Fig. [3] Errors in temperature are omitted for 
clarity. The dotted line indicates the solar value. 



In Fig. [3]four M 80 stars at Teif ~ 1 2 000 K stand out for their 
too high gravities. The fitting procedure outlined in ^3.11 was 
usually problematic for these targets, although both 600B and 
1400V spectra were of good quality and no problems were en- 
countered during reduction. Their spectra show a high quantity 
of strong metallic lines, but the fit with super-solar metallicity 
models is not convincingly better and, on the contrary, for two 
of these targets it was too poor to permit parameter determina- 
tion. For this reason, the tabulated and plotted results through- 
out this paper refer to fits with parent cluster metallicity mod- 
els, although they should not be the more appropriate. Their de- 
rived masses are very high compared to canonical values (2.5- 
5 Mq), and temperatures are completely inconsistent with their 
color: these stars can be found in Fig. [T] at the reddest end of 
the HB (U - V ~0.7), far from other HB stars with the same 
temperature and redder than the coolest "normal" HB stars with 
Teff =8 500 K. All these results point to a mismatch between their 
spectra and theoretical models, and we suspect these objects are 
foreground main sequence stars. This hypothesis is further cor- 
roborated by their metal-rich spectra, that would be very sur- 
prising for HB cluster members much redder than the Grundahl 
jump. Both temperatures and masses are roughly consistent with 
this explanation, but they might not be reliable. In particular, 
masses could be overestimated by Eq. [Tl due to a wrong dis- 
tance modulus. For one of these four stars we detect some RV 
variability on 1400V spectra, as analyzed in ^ 

Our results agree fairly well with theoretical expectations, 
and, within errors, the global behavior of all the points in Fig.[3]is 
to follow the canonical track, although the presence of some He- 
enhanced stars can not be excluded. One may note that the ten- 
dency toward lower gravities in the range Teff = 1 2 000- 1 7 000 K, 
that moves the points closer to the polluted models than the 
canonical ones, could also be explained by stellar winds unac- 
counted for by the model atmospheres. As discussed in M07, this 
hypothesis better explains the too low masses observed for many 
targets in this temperature range (see 34.11 ), because evolutionary 
effects or enhanced helium abundance would cause both lower 
gravity and higher luminosity, and the calculated mass would 
agree with expectations. 

4.1. Masses 

Our results o n ma sses resemble what found on NGC 6752 and 
discussed in IM07L and are plotted in Fig. |4] In brief, masses 



are systematically lower than theoretical prediction for stars 
cooler than 10000 K and in the range Teff=12 000-15 000 K, 
while they agree with models for Teff =15 000-23 000 K. Beyond 
Teff=23 000 K some stars are "normal", while others show an 
anomalously high mass. 

For temperatures below 10 000 K masses are too low in both 
clusters, as found in NGC 6752, despite the fact that the gravi- 
ties in Fig.[3]agree better with the canonical tracks. We conclude 
that the problem of low masses in this temperature range cannot 
be fully explained only by erroneously low gravities. Between 
15 000 and 23 000 K, masses fairly scatter around the theoreti- 
cal track. As already discussed in M07, this favors evolutionary 
effects or enhanced helium as explanation of the low gravities 
measured for many stars among 15 000 and 17 000 K, and it 
argues against the hypothesis of an enhanced stellar wind (un- 
accounted for in the models used for parameter measurements), 
that would cause underestimated masses. 

Above 23 000 K we cannot draw strong conclusions as done 
m [MOTI bee ause of the small number of stars and higher er- 
rors. In particular, in NGC 5986 we observed only five targets, 
and they are too scattered in the color-magnitude diagram to 
try to divide them in different families. On the other hand, in 
M 80 we can confirm what we found in NGC 6752. In fact, we 
overestimate masses with respect to theoretical predictions for 
four out of six stars, and all of them are systematically redder 
and slightly fainter than the two for which masses are "normal". 
These stars are indicated as full triangles in all figures. Hence, 
the strange dichotomy extensively discussed in M07 is present 
also in this cluster. We remind the reader that this behavior of 
masses is not a consequence of photometric data, because stars 
with higher derived masses are fainter, while a lower luminos- 
ity alone would imply a lower calculated mass. Hence, these 
stars appear both photometrically and spectroscopically distinct 
with respect to stars showing normal masses. The explanation 
for this behavior is still obscure. Higher masses are expected 
for EHB stars formed through a merging event of two He white 
dwarfs jHan et al J 12002 ^. but predictions are still much lower 
than measured values. We believe the too high masses are not 
physical, but just an effect of a mismatch between real stars and 
adopted model atmospheres. The same model atmospheres give 
good mass estimates for some stars and bad for others, suggest- 
ing that the actual stellar atmospheres are intrinsically different 
between the two groups of EHB stars. 
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Fig. 6. Maximum radial velocity variation observed for program stars. Thin errorbar indicates the 3cr interval. 



4.2. Helium abundance 

The heUum abundance was kept fixed to solar value for stars 
cooler than about 1 1 000 K, which showed no evidence of dif- 
fusion (see ^3.1b . Results for hotter stars are shown in Fig. |5] 
The plot reveals a clear trend with effective temperature fol- 
lowed by all star s in the three clusters, although it passed un- 
noticed by IMOTI because masked by the lack of stars between 
15 000 and 20000 K (see their Figure 8). The hehum abun- 
dance turns to sub-solar values at 12000 K, possibly reaching a 
minimum at 15 000-16000 K, then rises again steeply and con- 
tinuously up to about 22 000 K. Finally helium abundance de- 
creases again, but this transition region is undersampled. For 
stars hotter than 25 000 K the helium abundances scatter be- 
tween -3 < log ^ < -2. They possibly follow a double- 
peaked distribution rather than a wide single one, but this can be 
only guessed from the data. The same general trend can be seen 
also in Moehler et al. (2003, their Fig. 5) and, concerning the 
first decrease betwee n 12 000 and 15 000 K, in Behr (2003) and 
iFabbian e t al.' ('2005'), while it was probably hidden by larger er- 
rors in [Moehler et al. (2000). Helium depletion caused by diffu- 
sion is usually coupled with metal overabundances an d, interest- 
ingly e nough, even the iron abundances measured bv lPace et alj 
( l2006l their Fig. 4) in NGC 2808 follow the same trend observed 
here, reaching a maximum at about 15 000 K and then decreas- 
ing at higher temperatures. 

The observed behavior of helium surface abundances with 
temperature does not come unexpected. On the contrary, it con- 
firms theoretical expectations based on our current understand- 
ing of diffusion processes in the atmospheres of HB stars. Early 
results of Glaspev et al. (1989) and later detailed studies by 
many authors (see © showed that surface abundances abruptly 
change at Tgff ~11 000 K. The most common interpretation of 



the observed phenomena relies on the disappearance of the He i 
convection zone at this effective temperature, that would be re- 
sponsible for the gravitational settling of helium, with its con- 
sequent depletion in a thin radiative layer between the surface 
and the region of He ii ionization. Within this picture, this layer 
becomes thinner for increasing temper atures, because the second 
He ionization region moves outward ( Sweigart 2000^) , and more- 
over m ass-loss ( which contrasts diffus ion, Michaud & Charlandl 
il986h increases dVink & Cassisill2002[) . As a net result, the effi- 
ciency of diffusion would be expected to decrease with increas- 
ing temperature, as also suggested by the photometric properties 
of stars which, after the Grundahl jump, progressively reduce the 
discrepancies with standard models. This is exactly what we are 
observing between 15 000 and 23 000 K. Unfortunately, this sce- 
nario has recently been ruled out by the latest theoretical models, 
which show that diffusion affects layers much deepe r than the 
depth of the He i convection zone (IMichaud et al.l200 8, see their 
Figure 5). Actually, the models that best reproduce the observed 
trend of surface metal abundances assume that the outer regions 
down to j^ = 10"^ are completely mixed by turbulence, well be- 
low the He I ionization zone (Sweigart, priv. comm.). In other 
words, turbulence from this convecti ve layer is not able t o in- 
hibit levitation in Michaud's models. Michaud et al.l (l2008l) ten- 
tatively attribute the onset of levitation at a; 1 1 000 K to an abrupt 
change in stellar rotational velocity as observationally confirmed 
(e.g. Recio-Blanco et al.,.2002 : ,Behr.2003.) . but the issue is cur- 

ren tly under debate ; 

IMomanv et al.l ( l2002 l) proposed a new onset of diffusion at 
about 23 000 K, causing the photometrical anomaly known as 
"Momany Jump". Unfortunately the behavior of helium abun- 
dance among hotter stars is hard to decipher Both our present 
results and previous ones on NGC 6752 
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Table 4. Results for program stars. Columns 1: ID as in Table |2] Columns 2-7: absolute radial velocity, fundamental parameters 
(effective temperature, surface gravity, helium abundance and mass) and maximum radial velocity variation. 



ID 






logg log(l^) 


M 


A(RV)„,,, 




km s"' 


K 


Mo 


km s"' 


NGC 5986 



1 1 T 

1 IJVZ 


/ J±D 


1 /4UU±4UU 


4.30+U.U/ 


T 1 Q 1 n AA 


A /I 1 1 A AC 

U.41±U.Uj 


O 1 
O.l 


+ 


1 n 
i. 1 


17604 




O O CA 1 1 C A 

8s50±150 


T T C 1 A A"? 

3.25+0.07 


1 AA 1 A AA 

-1.00+0.00 


A y1 1 1 A AC 

U.41±U.Uj 








1 n^rt 1 




2 J /UO+1200 


5.5z+U.lz 


1 ii/i 1 n nn 

— i.o4+u.oy 


A '70 1 A 1 /I 

U. /o+U. 14 


111 

12.2 


+ 


5.1 


1 onnn 
loU/ / 


O 1 1 

oi±o 


y3jO±130 


3.38+0.14 


1 nn 1 n nn 
— l.UU+0.00 


r\ Al \ r\ rxn 
U.43+U.U/ 


A 1 

4. 1 


+ 


3. J 


1 /in 


/D+0 


o^cn 1 nn 
56j0±y0 


Q 1 n nc 
3.ZO+0.0J 


1 nn 1 n nn 
-l.UU+0.00 


A /I "7 1 A AC 


A 

4.0 


+ 


3. J 


ZlUi 


o/±j 


1 o 1 cn 1 1 zc\ 


Q OA 1 n nc 
3.04+U.0j 


— z.Zi+U.zl 


A /I C 1 A AC 


c o 

J.O 


+ 


3.Z 




on 1 Q 


1 nt^t\r\ 1 ^nn 
1 /dUU±4U0 


A A 1 I r\ r\n 
4.41+U.O / 


1 TO 1 n m 
— 1. /S+U.U/ 


A C/1 1 A AT 


1 1 A 

12. U 


+ 


A 1 

4. 1 




yZ± 1 


1 Tonn 1 Qnn 
1 /zU0±300 


A 1 n nc 
4.5j+U.0j 


1 m 1 n m 

— i.y /+U.U/ 


A C/1 1 A A/i 


o.y 


+ 


3. / 




0"^+ 1 1 
y 1 i 






Z. J JICU. i u 








7 n 


4930 


98±5 


10100+160 


3.48+0.07 


-1.00+0.00 


0.40+0.05 


4.8 


+ 


3.1 


5558 


84±6 


12500+150 


3.90+0.05 


-1.77+0.12 


0.60+0.07 


5.7 


+ 


3.3 


6102 


93±9 


15600+300 


4.30+0.05 


-2.37+0.14 


0.72+0.08 


8.6 


+ 


5.6 


7008 


81±5 


9000+160 


3.15+0.08 


-1.00+0.00 


0.32+0.04 


4.5 


+ 


3.2 


7430 


90±7 


9550+140 


3.37+0.07 


-1.00+0.00 


0.40+0.05 


9.1 


+ 


3.4 


8131 


89±7 


28300+700 


5.37+0.08 


-3.19+0.21 


0.62+0.08 


9.1 


+ 


6.0 


9049 


92±6 


18500+500 


4.67+0.07 


-1.70+0.09 


0.72+0.09 


13.1 


+ 


6.2 


9250 


91±6 


9150+140 


3.21+0.08 


-1.00+0.00 


0.33+0.04 


4.3 


+ 


3.3 


10390 


89±7 


9100+130 


3.12+0.07 


-1.00+0.00 


0.34+0.04 


5.4 


+ 


3.2 


11215 


91+7 


11000+110 


3.67+0.03 


-1.00+0.00 


0.51+0.05 


5.3 


+ 


3.3 


11571 


102±10 


27800+900 


5.53+0.12 


-2.71+0.21 


0.48+0.08 


20.4 


+ 


7.2 


12099 


103±11 


23700+800 


5.14+0.08 


-2.00+0.09 


0.52+0.07 


14.1 


+ 


5.9 


M80 


14327 


15±7 


20100+700 


5.00+0.08 


-1.44+0.05 


0.45+0.06 


8.1 


+ 


5.8 


14786 


6±4 


8280+50 


3.12+0.02 


-1.00+0.00 


0.45+0.05 


5.5 


+ 


3.3 


14985 


6±9 


17500+400 


4.78+0.05 


-2.25+0.09 


0.72+0.08 


12.6 


+ 


4.6 


15200 


12±6 


17800+500 


4.67+0.07 


-1.84+0.07 


0.60+0.08 


8.7 


+ 


4.3 


16389 


3±9 


18100+400 


4.75+0.05 


-2.00+0.07 


0.49+0.06 


23.1 


+ 


4.6 


16163 


5+5 


15400+200 


4.25+0.05 


-2.71+0.14 


0.46+0.05 


9.5 


+ 


4.6 


17173 


10±13 


12700+300 


4.23+0.12 


-1.00+0.00 


2.40+0.39 


26.4 


+ 


3.3 


17114 


8±4 


16500+300 


4.45+0.05 


-2.97+0.10 


0.78+0.09 


3.6 


+ 


3.3 


16707 


4±7 










11.1 


+ 


6.8 


17737 


1±5 


28100+700 


5.69+0.08 


-2.06+0.12 


0.91+0.12 


11.8 


+ 


6.1 


18516 


10±5 


18700+500 


4.78+0.07 


-2.00+0.07 


0.46+0.06 


8.2 


+ 


4.2 


18110 


10±4 


12400+150 


3.84+0.03 


-2.04+0.16 


0.48+0.05 


4.2 


+ 


3.1 


18992 


7±4 


8850+180 


3.16+0.08 


-1.00+0.00 


0.36+0.05 


2.9 


+ 


3.1 


19040 


15±5 


18800+400 


4.81+0.05 


-1.91+0.07 


0.51+0.06 


15.7 


+ 


6.2 


18391 


12±3 


12100+200 


4.37+0.07 


-1.00+0.00 


5.13+0.64 


5.7 


+ 


3.2 


19191 


13±9 


20600+600 


5.08+0.07 


-1.65+0.07 


0.49+0.07 


11.2 


+ 


5.1 


19127 


13±11 


25200+700 


5.51+0.08 


-2.94+0.17 


0.66+0.09 


18.1 


+ 


7.7 


12304 


7±6 


27300+600 


5.65+0.08 


-1.98+0.10 


0.71+0.10 


17.9 


+ 


7.3 


14201 


16±5 


16700+400 


4.43+0.05 


-2.34+0.12 


0.51+0.06 


2.4 


+ 


3.2 


13787 


1±4 


17700+300 


4.72+0.05 


-2.05+0.07 


0.63+0.07 


2.1 


+ 


3.2 


12663 


19±5 


27600+700 


5.38+0.08 


-3.14+0.17 


0.56+0.08 


16.5 


+ 


6.0 


12767 


20±9 


29900+500 


5.73+0.05 


-1.99+0.07 


0.66+0.08 


4.5 


+ 


4.7 


13839 


-4±11 


8450+40 


3.26+0.02 


-1.00+0.00 


0.55+0.06 


6.9 


+ 


3.2 


14022 


-11±9 


19900+800 


5.15+0.08 


-1.75+0.09 


0.76+0.11 


4.8 


+ 


4.9 


14387 


1±8 


15100+300 


4.36+0.05 


-2.37+0.14 


0.62+0.07 


7.3 


+ 


3.2 


12526 


6±15 


28000+700 


5.51+0.08 


-2.24+0.10 


0.52+0.07 


12.0 


+ 


6.3 


13179 


15±7 


16500+300 


4.17+0.05 


-1.34+0.05 


0.47+0.05 


3.5 


+ 


3.1 


15682 


14±10 


12400+200 


3.78+0.05 


-2.09+0.17 


0.45+0.05 


4.0 


+ 


3.1 


14327 


15±7 


20500+500 


5.06+0.07 


-1.42+0.05 


0.49+0.06 


5.8 


+ 


4.5 


15183 


0±15 


12600+400 


4.19+0.12 


-1.00+0.00 


3.09+0.06 


7.4 


+ 


3.2 


14813 


15±9 


8700+100 


3.16+0.05 


-1.00+0.00 


0.37+0.04 


3.1 


+ 


3.1 


15470 


8±6 


12600+300 


4.23+0.12 


-1.00+0.00 


2.95+0.48 


4.7 


+ 


3.2 
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dMoni Bidin et al.ll2007l: iMoehler et ani2000l) show that heUum 
is depleted between a factor of 10 and 100 for these stars, con- 
firming that diffusion is active at these temperatures, but with a 
scatter that is much larger than observational errors. Stars with 
the same temperature and different helium depletion seem to co- 
exist. Even evolutionary effects cannot be ruled out, as helium 
surface abundance could be a time-dependent result of compet- 
ing processes (Michaudet al. 1983). As already noted in M07, 
the helium abundances of EHB stars are not related with anoma- 
lous masses discussed in 14.11 For example, stars #12663 and 
#19127 are strongly depleted in hehum (log(N(He)/N(H))* -3), 
but the mass is "normal" for the first and too high for the second 
one. 



Th e comparison of our results with Figure 1 of lO'Toold 
(12008') can be very instructive. In that plot the author summa- 
rizes current knowledge about the trend of helium abundance 
with temperature for field sdB stars, gathering t h e data from 
many extensive survey s ( Edelmann et al.l 120031 : iLisker et"ar 



many extensive survey s ( Edelmann et al.l 120031 : iLisker et aL 
l2005tlStroer et al.l2005tlHirsch et al.ll2008h . The aim of lO' Toole 
(EoosT is to analyze the two families of field EHB stars dis- 



covered by Edelmann et al. (2003), one being a factor of ten 
more depleted in helium than the other Field stars cooler than 
20 000 K are undersampled, but they clearly do not follow the 
trend observed by us in three globular clusters. On the contrary, 
the helium abundances in the figure of O' Toole seem to decrease 
monotonically from 12000 K to at least 23 000 K. We have no 
explanation for this difference, but we think it deserves further 
investigation. 

At a first glance, results for globular cluster EHB stars agree 
fairly well with those for field EHB stars, with a possible small 
decrease of helium abundance between 20 000 and 23 000 K, and 
hotter stars being scattered between log(^) — -1.8 and -3.2. 
Hence, quite surprisingly, we would not be observing any evi- 
dence of a second, more depleted family of stars. These are only 
about 15% of field EHBs in the studied temperature range, and it 
might be possible to explain their absence in our sample by their 
relative scarcity alone. Nevertheless, we would statistically ex- 
pect about 4 such stars in our sample. If these He-poorer objects 
were stars with no core He-burning, as proposed by O' Toole 
(|2008|) . there would be no clear reason for their absence in globu- 
lar clusters. However, at a further inspection an alternative inter- 
pretation is possible. In fact, comparing the diagrams one could 
suspect that an offset is present for EHB stars, with our abun- 
dances being higher by about 0.5 dex. For example, our stars 
with Teff=2 0000 K cluster at about logfi^ = -1.5, while in 
Figure 1 of O'Toojd (120081) they are at about -2. Lowering our 
abundances by 0.5 dex in Fig.|5] the bulk of EHB stars would lie 
between log ^ = -2 and -3, while five stars would be much 
more depleted, at about log ^ - -3.7. In this case, the helium 
abundance of both the main population and He-depleted stars, 
and even the number ratio of the two families of EHB stars, 
would excellently agree with results among sdB stars. However, 
arguing for the presence of this offset is quite speculative, be- 
cause stars cooler than 20000 K show no evidence of it. The 
procedure adopted here to measure helium abundance is very 
standard for sdB studies, and we share the fitting routine with 
almost all fiel d surveys. Even the m odel atmospheres are often 
the same (e.g. lEdelmann et al.l2006h . and no offset has ever been 
observed. To verify an intrisic higher helium abundance for our 
cluster EHB stars would require a more extensive sample. As 
a conclusion, the presence of a helium-poorer EHB population 
in globular clusters, analogous to what is observed among field 
stars, must currently remain an open issue. 



Comparing our results with theoretical expectations is diffi- 
cult, because the effects of diffusion processes still lack a full 
comprehension. Recent calculations of Michaud et al. (2008!) 
well reproduce observed surface abundances, but their trend 
with effective temperature (that is our observational result) is 
still unexplored. Moreover, a detailed model prediction can- 
not neglect the counter- acting effects of stellar wind, which 
are still poorly known. In fact, they are suspected to play an 
important role, because helium abundances among sdB stars 
are much too high t o be accounted for by diffusion models 
(iMichaud et aL 19891). and weak stellar winds can explain th e 
discrepancv (" Fontaine & Chaveij|l99"7H Unglaub & Bues'' l998h . 
The combination of diffusion and stellar wind could also pro- 
duce time-dependent surface abundances, although on a rela- 
tively short time scale, that could affect the results by introducing 
an intrinsic star-to-star scatter. A detailed comparison of our re- 
sults with field sdB stars would also require a precise knowledge 
of the effects of metallicity. In fact, field sdBs should on aver- 
age be younger and more metal-rich than our sample. lUnglaubl 
(|2008|) showed that strong coupled stellar winds, involving H 
and He in addition to accelerated metals, are prevented for metal- 
poor EHB stars. The effects of this result on the helium surface 
abundance is not stated by the author, but it could be an impor- 
tant cl ue to interpret th e differences between our results and the 
plot of lO'Toold (l2008h . 

5. Radial velocity variations 

The arc lamp spectra for star #17604 in NGC5986 were dam- 
aged by lines of hot pixels, resulting in a wavelength calibration 
without the precision required by our aims. The target is very 
cool, adding little information to our search for binaries focused 
on EHB stars, and we considered safer to exclude it from analy- 
sis. 

Our radial velocity results are summarized in Fig. |6] where 
we plot the maximum variation for each star. Variations are taken 
always positive, the sign just being a consequence of the arbi- 
trary definition of "template" and "object" spectrum. We indi- 
cate the 3cr interval with thinner errorbars. 

We detect RV variations above 3cr for one of the four cool 
targets suspected to be foreground main sequence stars in M 80 
(®, indicating it could be a spectroscopic binary. We did not 
find evidence for a companion in our spectra, although at our 
low resolution its features could be easily hidden by the lines of 
the brighter primary. We found that all weak lines, in particu- 
lar Fe lines and the Mgib triplet, follow Hp in its RV variations. 
Therefore, the companion could be a compact object or a low- 
mass main sequence star. No other cool stars shows any variabil- 
ity within our errors. Hence, no close binary is detected among 
them, in agreement with recent results from a large sample of 
coo l HB stars in NGC6752 (104 HB stars with Teft <20000 K, 
Mo ni Bidin et al. 2008a). 

Star #14327 shows no noticeable RV variation in both field 
M 80a and M 80b. We also checked that no variation occurs be- 
tween its spectra in the two fields. Nevertheless, it is the only 
star observed seven times, and in order to avoid in the statistical 
analysis one star with different temporal sampling with respect 
to all the others, we will simply exclude results from slit 12b, 
thus considering this star as a normal target observed in field 
M 80a only. The choice of the slit to exclude is dictated by the 
better temporal sampling of M 80a field. 

Among the EHB stars we detect one interesting close 
binary candidate for each cluster, namely star #4175 in 
NGC5986 (maximum variation 26.9±7.0 km s"') and #16389 
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Table 5. Statistical analysis of RV variations for our target EHB 
stars. 



ID 


A(RV)n,,, 


cr 




p(< A(RV)„,J 




(km s-') 




(%) 


(%) 


NGC 5986 


17691 


12.2 


1.51 


73.6 


11.9 


4175 


26.9 


3.84 


0.15 


28.3 


8131 


9.1 


1.52 


73.6 


8.3 


11571 


20.4 


2.83 


4.0 


21.2 


12099 


14.1 


2.39 


13.4 


14.1 


M80 


14327 


8.1 


1.40 


60.4 


11.9 


14985* 


12.6 


2.74 


3.0 


19.5 


16389 


23.1 


5.02 


2.4- 10^' 


35.1 


16707 


11.1 


1.63 


44.6 


35.1 


17737 


11.8 


1.93 


26.7 


18.2 


19191 


11.2 


2.20 


14.2 


17.2 


19127 


18.1 


2.35 


10.0 


28.0 


12304 


17.9 


2.45 


3.5 


51.0 


12663 


16.5 


2.75 


1.5 


48.3 


12767 


4.5 


0.96 


65.2 


17.5 


14022 


4.8 


0.98 


65.2 


18.5 


12526 


12.0 


1.90 


14.4 


38.6 



in M80 (maximum variation 23.1±4.6 km s '). Variations for 
these targets are small compared to thos e observed on typ- 
ical sdB close binaries (s ee for example iMaxted et al.l 120011; 
iMorales-Rueda et aLll2003h . but so are the errors. In particular, 
CC errors for these measurements are very small compared with 
stars with similar S/N. In our experience, this could be an in- 
direct indication that the variation is real and not due to noise- 
induced distortion of the CC function. Some stars show varia- 
tions at the edge of the 3cr interval, although not formally outside 
it, and we want to verify if other candidates could hide among 
them. We performed statistical calculations on all EHB targets 
to translate these considerations into numbers. Results are sum- 
marized in Table |5] The star IDs and the maximum observed 
RV variations are indicated in column 1 and 2, respectively. In 
column 3 the maximum RV variation is given in units of cr. To 
estimate the significance of variations near the 3cr threshold, in 
column 4 we calculated the probability that a variation greater 
or equal to occurs among the measurements, assuming a 

normal distribution of errors. From basic statistical relations, we 
have: 



rA(RV)„ 



^ _ r f/A(RVWx|« 



(2) 



where n is the number of independent measurements (10 in 
NGC 5986, 6 in M 80a, and 3 in M 80b), and 



eTf(x) 



— e ' dt. 
\7r Jo 



(3) 



Finally, we calculated the fraction of typical sdB close binaries 
with random phase and inclination angle, a 0.5 companion, 
and period P<10 days, that in our survey would show a RV varia- 
tion not greater than the maximum observed. This is an estimate 
of how likely the small observed variations can really indicate a 
close binary. Results are tabulated in column 5. We included in 
the analysis one cool target (#14985 in M 80) with RV variations 
near 3cr, marked with an asterisk in Table|5] 

The numbers in Table |5] help to clarify the situation. First, 
the two stars that show variations greater than 3cr are confirmed 



to be good candidates, in view of the negligible probability of 
the observed datum being a random variation. Their variations 
are indeed small, but not enough to rule out binarity: about two 
out of three typical sdB binaries would show higher variations in 
our observations but, for example, this fraction would decrease 
if we assume a lower mass companion. 

For all other EHB targets, except for field M 80b because of 
poor temporal sampling, the variations are unlikely to be caused 
by binarity because of low probabilities in column 5, although it 
cannot be excluded on the basis of this datum alone. The num- 
bers in column 4 are more conclusive, and they do not allow us 
to claim the detection of any other binary candidate. In some 
cases the probability of the datum is low, but not negligible. One 
doubtful exception could be star #12663 in M 80, but the proba- 
bility of its observed RV variation being random is still at least 
one order of magnitude higher than for the two proposed can- 
didates. Moreover, it is important to notice that this probabil- 
ity is extremely sensitive to the adopted error value. We calcu- 
lated that probabilities lower than 4% in column 4 would double 
if errors were 10% larger. Despite our efforts, we feel that our 
error estimates are unlikely to be more precise than this value. 
Therefore these low probabilities should be considered accurate 
only to within a factor of two. 

In summary, the numbers support the hypothesis of binarity 
for only one star per cluster, the only ones which show variations 
greater than 3cr. This indicates that requiring a variation greater 
than 3cr for the detection of a binary system is a good choice, and 
in the statistical analysis we will assume this value as detection 
threshold. 

The only ot her EHB close binary di scovered to date in a 
globular cluster dMoni Bidin et al.l l2008al) sho ws a Mgib triplet 
anomalously strong for stars at its temperature dMoni Bidin et alj 
2007), indicative of a cool main sequence companion. We in- 
spected the sum of all the V1400 spectra, but we never found 
evidence for this feature in any EHB target. Therefore, the com- 
panions of our binary candidates are most likely compact objects 
(as white dwarfs), or very low-mass main sequence stars. In fact, 
we have no indication about the period, and very low-mass com- 
panions in close nearly edge-on orbits can c ause large RV varia - 
tions (see for example the NY Vir system, Vuck ovic et alj2007h . 
Thus, a low-mass main sequence companion is plausible. 

6. EHB close binary fraction 

In the statistical analysis of our results we will assume the de- 
tection of two EHB close binaries, one per cluster, as discussed 
in ^ The candidate in M 80 is slightly cooler than the typical 
temperature boundary for EHB stars (Tgff >20000 K). The dis- 
tinction is more than just a conventional definition, because it ap- 
proximately sets the transition between classical HB stars with 
envelopes sufficiently massive to sustain the hydrogen-burning 
shell, and EHB stars, which do not have this second energy 
source. The post-HB evolution is also very different, as depicted 
in fl] It can also be noted from Fig. [T] that at approximately 
20 000 K there is an underpopulated region in the HB of both 
clusters, and the binary in M80 is brighter than this gap (at 
V^18.7), indicating that it is indeed cooler than typical EHB 
stars. Despite these arguments, we consider that excluding the 
candidate from statistical analysis on the basis of its tempera- 
ture alone would be quite artificial. In fact, field surveys usually 
focus on hotter stars, excluding targets with lower temperatures 
(see for example Maxted et al. 2001), but they are sometimes 
included (Ulla & Theill 1998), and bina ry sys tems are detected 
among them (Aznar Cuadrado & Jeffervll2001h . Population syn- 
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Period (days) Binary fraction f 



Fig. 7. Probability of binary detection in our survey, as a function 
of binary period, solid line: NGC5986. dotted line: field M80a. 
dashed line: field M 80b. The calculation assumes a companion 
of 0.5 Mq in circular orbit. 



Fig. 8. Curves of probability for the close binary fraction /, as 
calculated from our results, solid line: NGC5986 (binaries with 
periods P<10 days), dotted line: M 80 (periods P<5 days). 



thesis models also indicate that sdBs as cool as 15 00 0-16000 K 
can b e formed by interactions within binary systems (iHan et al.l 
l2003h . Moreover, relatively massive (0.75 Mq) EHB stars, or 
with a relatively massive envelope (>0.01 Mq), move to temper- 
atures low er than 20 000 K during the first stages of post-EHB 
evolution dHan et al.ll2002l) . Hence, the candidate could even be 
an evolving object. Its exclusion would strengthen the results in 
favor of a lack of close systems. 

In the previous section we found that some hot stars show 
marginally significant RV variations. We emphasize that the fol- 
lowing analysis does not depend on whether we consider them as 
candidates or not. In fact, once the detection threshold is fixed, 
the only input is the number of stars with variations above and 
below it, and the routines automatically consider undetected sys- 
tems (because of unfavorable temporal sampling or low inclina- 
tion angles) as part of the calculations. 

The close binary detection probability of our survey was cal- 
culated as in Pa per L In brief, 2500 typical sdB binaries (with 
a 0.5 Mq companion in circular orbit, as assumed for example 
by iMaxted et all 120011: [Morales-Rueda et al.l l2006>) were simu- 
lated for each value of period P, evenly distributed in the phase- 
sin / space (where i is the inclination of orbit with respect to the 
line of sight). Then we calculated the fraction of these synthetic 
systems that would have been detected in our observations, i.e. 
showing RV variation greater than the detection threshold. We 
fixed the threshold at the 3cr value, because we found (© that 
variations lower than this value are not sufficiently significant. 
We adopted 3cr=20 km s ' for NGC5986 and 18 km s ' for 
M 80, which are average values for our program stars and well 
represent the typical accuracy of measurements. Results are plot- 
ted in Fig. |2] The detection probability for NGC 5986 is high, 
as a consequence of the good temporal sampling. On the con- 
trary, the lack of data in the first two nights for M 80, due to bad 



weather, strongly damaged the efficiency of the survey, mainly 
in the field M 80b. We will limit our analy sis t o perio ds P<5 days 
on this cluster, as done by lMoni Bidin et al.i (l2008al) for similar 
reasons. This limitation is not too severe, because our main aim 
is the comparison with results on field sdB stars. The exact shape 
of their period distribution is still unknown, but is suspected to 
be strongly peaked at P=lday, and binaries with P>5 days rep- 
resent only the tail of the d istribution (see for example Fig. 2 of 
iMorales-Rueda et al.ll2006l) . 

The probability of detecting Ns binaries out of a sample of 
N targets is: 

p = ^^^^ (dffm-Aff-^\ (4) 

^ (N-Nb)!Nb! ^ ' 

where / is the binary fraction and d is the probability of de- 
tection weighted with the period distribution. For M 80 we used 
the mean of the detection probability of the two fields, weighted 
with the number of EHB targets observed in each. The shape 
of the period distributio n affects the results only marginally, as 
akeady demonstrated by iPaper^ and [Moni Bidin et al. (2008ai). 
A Gaussian distri bution in log P cent ered on logP^O days, a s 
proposed by Maxted et al.l (1200 lb and iNapiwotzki et al.1 (|2004|) . 
does not change p(/) by more than 0.01-0.03 with respect to the 
flat distribution we assumed here. We used the relation (|4]l as a 
function of / to calculate the probability p(/), given our obser- 
vational results of one detection out of five targets in NGC 5986 
and one out of eleven in M 80. As stated before, we must limit 
our analysis to periods P<5 days on M 80. The results of our 
calculations are plotted in Fig. |8] The curves are far from being 
Gaussian-shaped, hence we cannot simply deduce a best esti- 
mate with an associated error Nevertheless, from their analysis 
we can draw important conclusions. 
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6.1. NGC5986 

In NGC5986 the most probable value of /p<iodays is 25%, as 
indicated by the peak of the curve in Fig. [8] This must be con- 
sidered the best estimate for the EHB close binary fraction in 
this cluster. Nevertheless, the small number of targets implies 
a very wide curve, so that no value of / can be safely ex- 
cluded . A very high binary fraction as proposed bv lMaxted et al.l 
(1200 1^ for field stars is not probable, but cannot be completely 
ruled out (p(/=0.7)=10.9%). The probability of our best esti- 
mate is p(/=0.25)=41%, which does not noticeably differ from 
p(/=0.4). Hence, a binary fraction as high as the lowest deter- 
mination among field sdB stars is not preferred, but it is per- 
fectly reasonable. There is no improvement in limiting to pe- 
riods shorter than 5 days, because the width of the probabil- 
ity curve is dictated by small number statistic and not by the 
sensitiveness of the survey. We conclude that our data do not 
give strong constraints on the EHB close binary fraction for 
this cluster, and they agre e both with the unexpected low esti- 
mates found in NGC 6752 (i Moni Bidin et al J2006bl,l2008a^ and 
with the lowest determinations for field sdBs dNapiwotz ki et aTl 
l20()4j). They tend to exclude high binary fractions as proposed 
bv lMaxted et al. (2001) for field sdBs. Given the most probable 
value, the results hint that the real fraction could be relatively 
low. 



6.2. M80 

Stronger results can be obtained on M80. The best estimate 
is /p<5days=12%. Very low compared to any determination 
among field stars. High values are ruled out (p(/=0.7)=0.2%). 
Re sults well agree with the extremely low fraction found 
bv 'Mon i Bidin et all (|2008a) in NGC 6752 (p(/=0.04)=25%), 
much more than the lowest values for field sdBs, which are very 
improbable for this cluster (p(/=0.4)=8.2%). Within a 90% con- 
fidence level, / is lower than 38%. 

In summary, although our results are not as strong as pre- 
vious ones on NGC 6752, we find that also in M80 EHB close 
binary systems are lacking, at variance with what is observed 
among field sdBs. The binary fraction is not very well con- 
strained, but values observed among field samples are very 
unlikely, w hile results agree well w ith the tiny 4% found in 
NGC 6752 dMoni Bidin et al.ll2008al) . M 80 is the second globu- 
lar cluster for which a lack of EHB close systems is found, and 
this indicates that it sh ould not be a pec uliarity of NGC 6752. 
Preliminary results by iMoni Bidin et al. I (f2008b) suggest that 
some clusters could be different, but the early stage of their anal- 
ysis and their small number statistics strongly call for further in- 
vestigation. In we discussed the models for the still unclear 
sdB star formation mechanisms. In light of these results, any 
successful model must take into account the significant differ- 
ence between stars inside and outside globular clusters. 



6.3. Comparison with Han (2008) results 



As mentioned in ^ lHanl(l2008h confirmed with theoretical cal- 
culations that the binar y scenario natural l y impli es a /-age re- 
lation, as proposed by iMoni Bidin et al.l (l2008ah . In brief, the 
model assumes that dynamical interactions in binary systems are 
responsible for sdB star formation, but the efficiency of the var- 
ious channels varies with the age of sdB progenitors, leading to 
a decreasing fraction of close binaries with increasing mean age 
of the population. 



The general conclusion of our work, that / is low also in a 
second globular cluster and possibly in a third one, is a fair con- 
firmation of Han's model. The fraction predicted in NGC 6752 
by his preferred set of parameters (2%) fairly agrees with ob- 
servational results o f|M oni Bidin et al.l(l2008al) . Unfortunately, a 
direct comparison for the two new globular clusters is compli- 
cated by the large uncertainties in both model predictions and 
empirical determinations. Therefore, the analysis cannot rely on 
exact values, but on the general behavior of the results. 

The best estimates of / in the three clusters so far surveyed 
are all smaller than those found for the field sdB stars, but too 
different from each other if compared to Han's prediction. The 
lowest observed fraction is for NGC 6752 and the highest is for 
NGC 5986. The order is coiTect, because NGC 6 752 is the oldest 
cluster while NGC 5986 is the youngest one dPe Angeli et alj 
2005i) . but age difference is too small, about 2 Gyr, compared 
with the modeled /-age relation. 

The parameter set preferred by Han creates an /-age relation 
that fairly agrees with observations of NGC 6752, but its predic- 
tions are not easily compatible with our new results: the model 
implies a nearly constant /=2% for populations older than 10 
Gyrs. Despite the uncertainties in our study, it is clear that this 
low value is improbable in both NGC 5986 (6% probability) and 
M 80 (14%, see Fig.[8]l. The combined probability, i.e. the prob- 
ability that / ^2% in both cluster is negligible, according to our 
observations. This tends to exclude the /-age relation derived by 
this model. Other models studied by Han (2008, see his Fig. 3), 
with different sets of input parameters, predict higher close bi- 
nary fractions and a steeper relation with age, thus solving the 
outlined contradictions, but their high expected / are incompat- 
ible with measurements in NGC 6752, the more robust of the 
observed results. 

From this analysis we conclude that, despite the good general 
agreement (close EHB binaries are predicted and observed to be 
a minor population), models and observations still lack a good 
agreement on the details, although observational constraints are 
still not strong enough to be conclusive. Maybe the discrepancies 
could be mitigated by a refined set of model parameters, or some 
other effect (like dynamical interactions in dense environments) 
could be invoked to slightly change / from cluster to cluster. 

It is important to note that the present comparison strongly 
relies on our results on M 80, which are the most precise, and 
the low temperature of the binary candidate in this cluster leaves 
spac e for doubts on the feasibility of the comparison itself. In 
fact. iHanl (1200 8) did not apply any temperature cut, but his /-age 
relation was obtained taking into account a GK selection effect, 
and models corrected to con sider this bias tend to exclude sdBs 
as cool as our target dHan et al. 2003). Anyway, the system in 
M 80 should not be a EHB-hMS wide binary (© as the systems 
selected against by the GK effect, so the issue remains uncer- 
tain. Han (2008) just states that, in absence of the GK effect, his 
predicted / should be smaller at any age. 

We conclude that the confirmation and refinement of iHanl 
(I2008h model require more precise empirical / measurements in 
these and other globular clusters, that could even help constrain- 
ing the model parameters, in particular the physically-impo rtant 
and poorly-known common envelope efficiency acE (see iHanl 
i2008L for a discussion). 



' The GK selection effect is an observational bias against EHB stars 
with a companion of G-K spectral type, because of the composite spec- 
trum, or earlier, because of spectral dominance of the companion. 
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7. Summary 

We analyzed radial velocity variations for 5 1 hot HB/EHB stars 
in two Galactic globular clusters, in search for signatures of close 
binary systems. We also studied low-resolution spectra of pro- 
gram stars measuring temperature, surface gravity, helium abun- 
dance and mass. Our main results can be summarized as follows: 

- In M 80 we confirm the ano malous behavior of spec troscopic 
masses found in NGC 6752 (iMoni Bidin et al.l2007l) for stars 
hotter than 23 000 K, although this result is less evident 
due to a smaller sample. Stars being fainter and/or redder 
show too high masses with respect to theoretical expecta- 
tions, whereas their lower luminosities would suggest lower 
masses. In NGC 5986 the small number of EHB target ob- 
served prevents such an analysis. 

- For the first time we observe a clear trend of helium abun- 
dance with temperature along the entire blue HB. These 
results confirm that helium depletion due to atmospheric 
diffusion reaches a maximum at about 15 000 K, then the 
helium abundance rises again with increasing temperature. 
This behavior agrees qualitatively with expectations for the 
effects of diffusion at different temperatures along the HB. 
Somewhere at about 23 000 K the helium abundance could 
start to decrease again, but the observed pattern is hard to 
decipher. For these hot stars helium is depleted between a 
factor of 10 and 100, without a clear trend with temperature, 
nor a relation with the previously mentioned dichotomy on 
calculated masses. 

- We detect one EHB close binary candidate per cluster. Their 
RV variations are quite small compared to typical sdB bi- 
nary systems, but the probability of their being due to ran- 
dom errors is negligible. The candidate in M 80 is slightly 
cooler than typical EHB stars (Teft-18 100 K). None of 
them, nor any other EHB target, show the Mgib triplet as 
signature of a cool companion, at variance with the only 
other sdB close binary dis covered to date in a globular cluster 
dMoni Bidin et alj|2008ah . Therefore their companions are 
more likely compact objects such as white dwarfs, or very 
low-mass main sequence stars. 

- The best estimate for the EHB close binary fraction in 
NGC 5986 is /=25%. This suggests that the fraction could 
be small, but no value lower than 70% can safely be ex- 
cluded because of the small observed sample. Nevertheless, 
the probability of a very high binary fraction is low. 

- In M 80 EHB close binaries with period shorter than 5 days 
are lacking with respect to what is observed among field sdB 
stars. The best estimate is /p<5days = 12%, and within a 90% 
confidence level /p<5days <38%. The fraction is not very well 
determined, but even the lowest values found for field stars 
are very improbable and can be ruled out. R esults fairly agree 
with th e tiny /=4% found in NGC 6752 ( Moni Bidin et alj 
l2008ah . M80 is the second cluster for which this behavior 
is observed. This indicates that it is not just a peculiarity 
of NGC 6752 but it should be quite a common feature, pos- 
ing important new constraints on models concerning sdB star 
formation. 

- Our results agree with the existence of an /-age relation 
for sdB star populations, as proposed by Moni Bid in et al] 
(l2008ah and theoretically modeled by Han. (20081) in 
the framework of binary scenarios for their formation. 
Uncertainties in both predicted and observed values prevent 
their direct comparis on. N e verthe less, we find that the /-age 
relation modeled by iHanl (|2008|) has some problem in re- 
producing the details of the observations. A refined set of 



model parameters could be needed, or the intervention of 
some other effect able to affect / and varying from cluster 
to cluster, like for instance dynamical interactions in dense 
environment. 

The low binary fraction among globular cluster EHB stars is 
consistent with the typical (low) binary fraction among globu- 
lar cluster stars. After all, the high binary fraction of sdB stars 
might not be a peculiarity related to their nature of hot He burn- 
ing stars, but a simple consequence of the high binary fraction 
among field stars. If so, we still need to identify the cause of 
sdB star formation. Nevertheless, recent refined binary models 
can (at least qualitatively) account for the lack of close EHB 
systems in globular clusters, still retaining the hypothesis of a 
binary origin (Han 2008). Thus, the "binary scenario" cannot be 
ruled out by current observations, and we are still far from a full 
understanding of the complete picture. 
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